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Masses of tetraquarks with two heavy quarks and open charm and bottom are cal- 
culated in the framework of the diquark-antidiquark picture in the relativistic quark 
model. All model parameters were regarded as fixed by previous considerations of 
various properties of mesons and baryons. The light quarks and diquarks are treated 
completely relativistically. The c quark is assumed to be heavy enough to make the 
diquark configurations dominating. The diquarks are considered not to be point- 
like but to have an internal structure which is taken into account by the calculated 
diquark form factor entering the diquark-gluon interaction. It is found that all the 
{cc){qq') tetraquarks have masses above the thresholds for decays into open charm 
mesons. Only the I(J^) = 0(1+) state of {bb){ud) lies below the BB* threshold and 
is predicted to be narrow. 
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I. INTRODUCTION 



Recent experimental studies of the heavy meson spectroscopy revealed several new states, 
such as X(3872), F(4260), D!q(2317) etc., which cannot be simply accommodated in the 
quark- ant iquark (qq) picture [l[ . These states can be considered as indications of the possible 
existence of exotic multiquark states which were proposed long ago, e.g. in The idea 
to revisit the multiquark picture using diquarks has been put forward by Jaffe and Wilczek 
[3|. At present, vast experimental and theoretical evidence of the important role played by 
diquark correlations in hadrons is collected 
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The simplest multiquark system is a tetraquark, composed of two quarks and two anti- 
quarks. Heavy tetraquarks are of particular interest, since the presence of a heavy quark 
increases the binding energy of the bound system and, as a result, the possibility that such 
tetraquarks will have masses below the thresholds for decays to mesons with open heavy 
flavour. If such strong decays are kinematically forbidden, then corresponding tetraquarks 
can decay only weakly or electromagnetically and thus they should have a small decay 
width. In this paper we consider tetraquarks with two heavy quarks as bound systems of a 
diquark and antidiquark. Therefore we assume that both b and c quarks are heavy enough 
to make the attractive QQ interaction stronger than the Qq one. For tetraquarks contain- 
ing c quarks the obtained results crucially depend on this assumption. In particular, the 
doubly heavy {QQ'){qq') tetraquark {Q = b,c and q = u,d,s) is considered as the bound 
system of the heavy diquark {QQ') and light antidiquark {qq'), while the {cq){bq') tetraquark 
is the bound state of the heavy- light diquark (eg) and antidiquark {bq'). Masses of heavy 
tetraquarks with hidden charm {cq){cq) and bottom {bq){bq) were calculated in our previ- 
ous paper [5]. There the dynamical analysis has shown that X(3872) and y(4260) can be 
indeed the diquark-antidiquark tetraquarks with hidden charm. It was also argued that the 
corresponding ground-state tetraquarks with hidden bottom have masses below the open 
bottom threshold, and thus they should be narrow states. 

It is important to investigate the possible stability of the {QQ'){qq') tetraquarks since 
they are explicitly exotic states with the heavy flavour number equal to 2. Thus, their 
observation would be a direct proof of the existence of the multiquark states. Estimates of 
the production rates of such tetraquarks indicate that they could be produced and detected 
at present (SELEX, Tevatron, RHIC) and future facilities (LHC, LHCb, ALICE) Q. 

To calculate the masses of heavy tetraquarks we use the relativistic quark model based on 
the quasipotential approach in quantum field theory. Previously we considered in our model 
the mass spectra of the ground-state and excited doubly- heavy {QQq) [2] and heavy {qqQ) jsl, 
0] baryons in the heavy-diquark-light-quark and light-diquark-heavy-quark approximations, 
respectively. The light quarks and light diquarks were treated completely relativistically. 
The internal structure of the light and heavy diquarks was taken into account by calculating 
diquark-gluon form factors in terms of the obtained diquark wave functions. Such scheme 
proved to be very effective and successful in our calculation of the masses of heavy baryons 
in good agreement with experimental data [lo|. The predicted masses of the T^b, and 
Sfc baryons proved to be very close to the recently measured ones [HI, H Q. Moreover, 



in Ref. [Qj it was shown that currently available experimental data on excited charmed 
baryons can be accommodated in the picture treating a heavy baryon as the bound system 
of the light diquark and heavy quark, experiencing orbital and radial excitations. This gives 
us additional confidence in the reliability of the diquark approximation within our model 
and motivates the consideration of tetraquarks as diquark-antidiquark bound systems. It is 
important to note that all parameters of our model were fixed in the previous calculations 
of meson mass spectra and decays, and we will keep their values in the following analysis of 
heavy tetraquarks. 



II. RELATIVISTIC QUARK MODEL 



In the quasipotential approach and diquark-antidiquark picture of heavy tetraquarks 
the interaction of two quarks in a diquark and the diquark-antidiquark interaction in a 
tetraquark are described by the diquark wave function (\l/rf) of the bound quark-quark state 
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and by the tetraquark wave function of the bound diquark-antidiquark state, respec- 
tively, which satisfy the quasipotential equation of the Schrodinger type [3] 

where the relativistic reduced mass is 

_ E,E, _ - (ml - ml)^ 



and El, E2 are given by 



IVP -ml + ml M"^ -ml + ml 
E-\ = , Eo = . (3) 

Here, M = Ei + E2 the bound-state mass (diquark or tetraquark), mi 2 are the masses of 
quarks (gi and ^2) which form the diquark or of the diquark {d) and antiquark {d') which 
form the heavy tetraquark (T), and p is their relative momentum. In the center-of-mass 
system the relative momentum squared on mass shell reads 

2 ^ [M^ - (mi + m^YWM^ - (mi - m^f] 

The kernel l^(p, q; M) in Eq. ([1]) is the quasipotential operator of the quark-quark or 
diquark-antidiquark interaction. It is constructed with the help of the off-mass-shell scatter- 
ing amplitude, projected onto the positive-energy states. In the following analysis we closely 
follow the similar construction of the quark-antiquark interaction in mesons which were ex- 



tensively studied in our relativistic quark model [IJ, |15| . For the quark-quark interaction in 
a diquark we use the relation Vqq = Vgg/2 arising under the assumption of an octet structure 
of the interaction from the difference in the qq and qq colour states.^ An important role 
in this construction is played by the Lorentz structure of the confining interaction. In our 
analysis of mesons, while constructing the quasipotential of the quark-antiquark interaction, 
we assumed that the effective interaction is the sum of the usual one-gluon exchange term 
and a mixture of long-range vector and scalar linear confining potentials, where the vector 
confining potential contains the Pauli terms. We use the same conventions for the con- 
struction of the quark-quark and diquark-antidiquark interactions in the tetraquark. The 
quasipotential is then defined as follows 0, 15 



(a) For the quark-quark {qq'), (Qq), {QQ') interactions, V{p,q; M) reads 

V{p, q; M) = Ui{p)u2{~p)V{p, q; M)ui{q)u2{-q), (5) 

with 



V(p,q;M) = i 



^asD,Ml^l2 + Krnf(k)r^(k)r2;,(-k) + Kif(k) 



Here, is the QCD coupling constant; D^^, is the gluon propagator in the Coulomb gauge, 
D''{k) = -^, D^^{k) = -^(6^^-^], D°^ = D^° = 0, (6) 



^ Obviously, it is important to study diquark correlations in gauge-invariant color-singlet hadron states on 
the lattice. 
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and k = p — q; 7^ and u{p) are the Dirac matrices and spinors 

u\p) 



e{p) + m 



\ Hp) 



1 

cr ■ p 



X 



(7) 



e{p) + m 



with e{p) = Vp^ + m'^. 

The effective long-range vector vertex of the quark is defined 3] by 



r^(k) = 7^ + ^(^,.uk' 



k = {0,k), 



(8) 



where k is the Pauh interaction constant characterizing the anomalous chromomagnetic 
moment of quarks. In configuration space the vector and scalar confining potentials in the 
nonrelativistic limit reduce to 



(l-£)Konf(r), 



with 



Ar + B, 



confv / ' '^confl 

where e is the mixing coefficient. 

(b) For the diquark-antidiquark {dd!) interaction, l^(p, q; M) is given by 

^d^d <J ^yi^d'^d' 

+rAPWd'{n \Jd,,r,,v^oA^) + Kif(k)l MQ)MQ') 



(9) 
(10) 



where {d{P)\J fj\d{Q)) is the vertex of the diquark-gluon interaction which takes into account 
the finite size of the diquark and is discussed below P*^'^ = (-E^o, ±p) and Q^'^ = {E^o, ±q), 
Ed = (M2 - M| + MJ)/(2M) and Ed' = (M^ - M| + M|)/(2M) " . 

The diquark state in the confining part of the diquark-antidiquark quasipotential ffTTl) is 
described by the wave functions 

'1 for a scalar diquark, 

^ Edip) for an axial- vector diquark. 



(12) 



where the four-vector 

£d{p) = 



[Sd ■ P) 

Md 



,Sd + 



[^d ■ P)P 



Md{Ed{p) + Md) 



^d{p)Pi^ = 0, 



(13) 



is the polarization vector of the axial- vector diquark with momentum p, Ed{p) = y P^ + M^, 
and erf(O) = (0, e^) is the polarization vector in the diquark rest frame. The effective long- 
range vector vertex of the diquark can be presented in the form 



J, 



2\ EdEd 



d;ti 



^JEdEd 



2Md 



for a scalar diquark, 
for an axial-vector diquark. 



(14) 



5 



TABLE I: Masses M and form factor parameters (for definitions see Eq. (jl9p ) of light diquarks. 
S and A denote scalar and axial-vector diquarks, antisymmetric [q, q'] and symmetric {q, q'} in 
flavour, respectively. 



Quark content Diquark type M (MeV) ^ (GeV) C (GeV^) 

[u,d] S 710 1.09 0.185 

{n,4 A 909 1.185 0.365 

[u,s] S 948 1.23 0.225 

{m,s} a 1069 1.15 0.325 

{s,s} A 1203 1.13 0.280 



where k = (0,k). Here, the antisymmetric tensor T,'^ is defined by 

iJ:,^)l = -tig,,5:-g,j;), (15) 

and the axial-vector diquark spin is given by {Sd-k)ii = —iskiu l^d is the total chromo- 
magnetic moment of the axial-vector diquark. 

The constituent quark masses rrih = 4.88 GeV, rric = 1.55 GeV, = = 0.33 GeV, 
rus = 0.5 GeV and the parameters of the linear potential A = 0.18 GeV^ and B = —0.3 GeV 
have the values typical in quark models. The value of the mixing coefficient of vector 
and scalar confining potentials e = —1 has been determined from the consideration of 



charmonium radiative decays [IJ] and the heavy-quark expansion [16|. The universal Pauli 
interaction constant k = —1 has been fixed from the analysis of the fine splitting of heavy 
quarkonia ^Pj- states jl3]- In this case, the long-range chromomagnetic interaction of quarks 
vanishes in accordance with the flux-tube model. 



III. LIGHT, HEAVY-LIGHT AND HEAVY DIQUARKS 

As the first step, we calculate the masses and form factors of the diquarks. As is well 
known, the light quarks are highly relativistic, which makes the v/c expansion inapplicable 
and thus a completely relativistic treatment is required. To achieve this goal in describing 
light and heavy-light diquarks, we closely follow our recent consideration of the spectra of 
light mesons and adopt the same procedure to make the relativistic quark potential local 
by replacing ei_2(p) = \l^i^2 + by £'1^2 (see the discussion in Ref. The resulting 

quark-quark interaction potential is equal to 1/2 of the qq interaction in the meson. We 
solve numerically the quasipotential equation with this complete relativistic potential which 
depends on the diquark mass in a complicated, highly nonlinear way. The obtained ground- 
state masses of scalar and axial- vector light and heavy diquarks [HI, 0, 0] are presented in 
Tables [T] and [m 

In order to determine the diquark interaction with the gluon field, which takes into ac- 
count the diquark structure, it is necessary to calculate the corresponding matrix element 
of the quark current between diquark states. This (diagonal) matrix element can be param- 
eterized by the following set of elastic form factors: 

(a) scalar diquark (5) 



{S{P)\J,\S{Q)) = h+{k'){P + Q),, 



(16) 
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TABLE II: Masses M and form factor parameters (for the definitions see Eq. (jl9|) ) of heavy- light 
and heavy diquarks. S and A denote scalar and axial-vector diquarks, antisymmetric [Q, q] and 
symmetric {Q, q} in flavour, respectively. 



1 1 1 1 O 

ILdL K 

content 


\ It /"i 1 1 o y \^ 

Uii^lLdbL is. 

type 




o — 

— c 






D — h 




M (MeV) 


e (GeV) 


C (GeV2) 


M (MeV) 


e (GeV) 


C (GeV2) 


[Q,u] 


S 


1973 


2.55 


0.63 


5359 


6.10 


0.55 


{Q,u} 


A 


2036 


2.51 


0.45 


5381 


6.05 


0.35 


[Q,s] 


S 


2091 


2.15 


1.05 


5462 


5.70 


0.35 


{Q,s} 


A 


2158 


2.12 


0.99 


5482 


5.65 


0.27 


[Q,c] 


S 








6519 


1.50 


0.59 


{Q,c} 


A 


3226 


1.30 


0.42 


6526 


1.50 


0.59 


{Q,b} 


A 


6526 


1.50 


0.59 


9778 


1.30 


1.60 



(b) axial-vector diquark (A) 

{A{P)\J,\A{Q)) = -[el{P).e,{Q)]h{e){P + Q), 

+h,{e) {[em ■ Q]e,,,{Q) + [e,{Q) ■ P]4,(P)} 

+h3{k')^Js:iP) ■ Q][e,{Q) ■ P]{P + g)„ (17) 

where k = P — Q, and Sd{P) is the polarization vector of the axial- vector diquark (fT3|) . 

Using the quasipotential approach with the impulse approximation for the vertex function 
of the quark-gluon interaction, we find [1] 

hik') = 0, 



Ea + Mj (27r 



r^d P + 



2e2(p) 



eiip + k) + ei{p) 



P 



Ed + Md J \ ti{p + k)+mi [2^ei{p + k)ei{p) 

vl/rf(p) + (1 ^ 2), (18) 



F(k^) 



2yei(p + k)ei{p){ei{p) + mi). 

where are the diquark wave functions. We calculated the corresponding form factors 
F{r)/r, which are the Fourier transforms of F(k^)/k^, using the diquark wave functions 
found by numerically solving the quasipotential equation. Our estimates show that this 
form factor can be approximated with high accuracy by the expression 



F(r) = 1 - e 



(19) 



which agrees with previously used approximations [7|. The values of the parameters ^ and C. 
for light, heavy-light and heavy scalar diquark [q, q'] and axial- vector diquark {q, q'} ground 
states are given in Tables HI and HTl 
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IV. MASSES OF HEAVY TETRAQUARKS 

As the second step, we calculate the masses of heavy tetraquarks considered as bound 
states of diquark and antidiquark. For the potential of the diquark-antidiquark interaction 
(fTTD we eet ^ 



V{r) = Vboui(r) + Konf(r) + -< 



M2{E2 + M2) 



VL 



conf \ 




+ 



Ei{Ei + Mi) E2{E2 + M2 



+ 



2 
1 



+ 



1 \ K'o^nf(r) 



Ml M| 



>L- (S1 + S2) 



^2(^2 + M2 

1 \ V!^ 



Ml 




L- (S 



Ml (El + Ml 



Mi(Ei + Mi) M2{E2 + M2) 



E1E2 



P 



Vcouiir) + K 



conf \ 



1 

+ - 

r 



Ud 

4 



+- 



Ei_ 
Ml 



4 VMi 

E^- ^^'^ 
M2 



(Si - S2 



-v;^oui(^) 



V^c'oul(0 + 



^EiE2_ 
4 M1M2 



»+^Coul( 

(S1 + S2) 



r) 



2r 



confl'") ^onf( 



X 



+ - 



:Si-r)(S2-r)-Si-S2 



Si • S2 



(20) 



where 



4 Fi(r)F2(r) 
3 r 



is the Coulomb-like one-gluon exchange potential which takes into account the finite sizes 
of the diquark and antidiquark through corresponding form factors Fi 2(r). Here, Si 2 and 
L are the spin operators of diquark and antidiquark and the operator of the relative orbital 
angular momentum. Since we limit our considerations to the ground states of heavy diquark- 
antidiquark bound systems ((L^) = 0), the spin-orbit and tensor terms in the potential (120|) 
do not contribute in the further analysis. In the following we choose the total chromomag- 
netic moment of the axial- vector diquark = 0. Such a choice appears to be natural, since 



^ In our paper [5| first two spin-orbit terms were missed. However they do not influence published numerical 
results, since masses mostly of ground states were calculated. Orbital excitations were considered only for 
the tetraquarks composed of the scalar diquark and scalar antidiquark for which the missed terms vanish. 
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TABLE III: Masses M of heavy-diquark (QQO^light-antidiquark {qq) states. T is the lowest 
threshold for decays into two heavy-light (Qq) mesons and A = M—T. All values are given in MeV. 



System 


State 




Q = Q' = c 




Q 


= Q' = b 




Q 


= c, Q' 


= b 






M 


T 


A 


M 


T 


A 


M 


T 


A 


iQQ'){ud) 


0(0+) 














7239 


7144 


95 




0(1+) 


3935 


3871 


64 


10502 


10604 


-102 


7246 


7190 


56 




1(1+) 














7403 


7190 


213 




1(0+) 


4056 


3729 


327 


10648 


10558 


90 


7383 


7144 


239 




1(1+) 


4079 


3871 


208 


10657 


10604 


53 


7396 


7190 


206 




1(2+) 


4118 


4014 


104 


10673 


10650 


23 


7422 


7332 


90 


iQQ'){us) 
















7444 


7232 


212 






4143 


3975 


168 


10706 


10693 


13 


7451 


7277 


174 


















7555 


7277 


278 






4221 


3833 


388 


10802 


10649 


153 


7540 


7232 


308 






4239 


3975 


264 


10809 


10693 


116 


7552 


7277 


275 




^(2+) 


4271 


4119 


152 


10823 


10742 


81 


7572 


7420 


152 


iQQ')iss) 


0(1+) 














7684 


7381 


303 




0(0+) 


4359 


3936 


423 


10932 


10739 


193 


7673 


7336 


337 




0(1+) 


4375 


4080 


295 


10939 


10786 


153 


7683 


7381 


302 




0(2+) 


4402 


4224 


178 


10950 


10833 


117 


7701 


7525 


176 



the long-range chromomagnetic interaction of diquarks proportional to jjd then also vanishes 
in accordance with the flux-tube model. 

We substitute the quasipotential fl2U]) in the quasipotential equation ([T]) and solve the 
resulting differential equation numerically. The calculated masses M of tetraquarks with 
open charm and/or bottom composed from the heavy diquark, containing two heavy quarks 
{QQ', Q = b,c), and the light antidiquark {qq', q = u,d,s) are presented in Table UTTl In 
this table we give the values of the lowest thresholds T for decays into two corresponding 
heavy-light mesons [{Qq) = D^*\ Di*\ B^*\Bi*^], which were calculated using the measured 
masses of these mesons [lo[. We also show values of the difference of the tetraquark and 
threshold masses A = M — T. If this quantity is negative, then the tetraquark lies below 
the threshold of the decay into mesons with open flavour and thus should be a narrow state 
which can be detected experimentally. The states with small positive values of A could be 
also observed as resonances, since their decay rates will be suppressed by the phase space. 
All other states are expected to be very broad and thus unobservable. We find that the only 
tetraquark which lies considerably below threshold is the 0(1+) state of {bb){ud). All other 
{QQ'){qq') tetraquarks are predicted to lie either close to or significantly above corresponding 
thresholds. It is evident from the results presented in Table IIIII that the heavy tetraquarks 
have increasing chances to be below the open flavour threshold and thus have a narrow width 
with the increase of the ratio of the heavy diquark mass to the light antidiquark mass. 



9 



TABLE IV: Masses M of diquark (cg')-antidiquark [hq) states. T is the lowest threshold for decays 
into two heavy-light {Qq) mesons and A = M — T; T' is the threshold for decays into the Sc*^ and 



a light meson [q'q) 


and A' 


= M - 


T'. All values are 


given 


in MeV. 










System 


State 






q' = u 










q' = s 






M 


T 


A 


r 


A' 


M 


T 


A 


r 


A' 


icq'Mhu) 


























0+ 


7177 


7144 


33 


6818 


359 


7294 


7232 


62 


6768 


526 




1+ 


7198 


7190 


8 


6880 


318 


7317 


7277 


40 


6820 


497 




1+ 


7242 


7190 


52 


6880 


362 


7362 


7277 


85 


6820 


542 




0+ 


7221 


7144 


77 


6818 


403 


7343 


7232 


111 


6768 


575 




1+ 


7242 


7190 


52 


6880 


362 


7364 


7277 


87 


6820 


544 




2+ 


7288 


7332 


-44 


7125 


163 


7406 


7420 


-14 


7228 


178 


{cq'){bs) 


























0+ 


7282 


7247 


35 


6768 


514 


7398 


7336 


62 


6818 


580 




1+ 


7302 


7293 


9 


6820 


482 


7418 


7381 


37 


6880 


538 




1+ 


7346 


7293 


53 


6820 


526 


7465 


7381 


84 


6880 


585 




0+ 


7325 


7247 


78 


6768 


557 


7445 


7336 


109 


6818 


627 




1+ 


7345 


7293 


52 


6820 


525 


7465 


7381 


84 


6880 


585 




2+ 


7389 


7437 


-48 


7228 


161 


7506 


7525 


-19 


7352 


154 



In Table HVl the calculated masses M of tetraquarks composed from a (eg) diquark and a 
{bq) antidiquark are listed. ^ We also give the lowest thresholds T for decays into heavy-light 
mesons as well as thresholds T' for decays into the S^*-* and light (q'q) mesons and A'- '' = 
]\^ _ )_4 gj-^^ ^j^g^i- Qj^iy 2+ states of {cq')(bq) have negative values of A and thus they 
should be stable with respect to decays into heavy-light {B and D) mesons. The predicted 
masses of lowest l"*" states of {cu)(bu) and {cu)(bs) tetraquarks lie only slightly above the 
corresponding thresholds T. However, all {cq){bq) tetraquarks are found to be significantly 
above the thresholds T' for decays into the -B^*-* and light (q'q) mesons. Nevertheless, the 
wave function of the spatially extended (cqMbq) tetraquark would have little overlap with 
the wave function of the compact Be meson 1^, thus substantially suppressing the decay 



rate in this channel. Therefore the above-mentioned {cq){bq) tetraquark states which are 
below the BD threshold have good chances to be rather narrow and could be detected 
experimentally. 

It is important to note that the comparison of the masses of heavy tetraquarks given in 
Tables UTTl and ITVl with our previous predictions [5| shows that the {QQ'){qq') states are, 
in general, heavier than the corresponding {Qq){Q'q') ones. This result has the following 
explanation. Although the relation [13] Mqq + < 2MQq holds between diquark masses, 
the binding energy in the heavy-light diquark (Qq)- heavy-light antidiquark (Qq) bound 



^ Such (cm) (bu) tetraquarks were recently argued [18| to be the best candidates for experimental detection. 
^ For the non-strange {cq)(bq) tetraquarks we give thresholds T' for decays of the / = states into B^*^ 

and T] or lu. These states should be more stable than the I — 1 ones, since their decays to Be and tt 

violate isospin. 
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TABLE V: Comparison of different theoretical predictions for the masses of heavy {QQ'){qq') 
tetraquarks (in MeV). 



System 


I{JP) 


this work 


[21] 


[22] 


[23] 


m 


[25] 


[26] [27] 


(cc) (ud) 




















0(1+) 


3935 


3931 


3876 




3764 


3927 


3905 4000 ± 200 




1(0+) 


4056 








4150 








1(1+) 


4079 








4186 








1(2+) 


4118 








4211 






{bb){ud) 




















0(1+) 


10502 


10525 


10504 


10558 


10261 


10426 


10200 ± 300 




1(0+) 


10648 




10587 


10766 


10690 








1(1+) 


10657 


10712 


10644 


10774 


10698 








1(2+) 


10673 


10735 




10790 


10707 






{bb)ius) 




10706 
10823 


10680 
10816 












{bc){ud) 


0(0+) 
0(1+) 
1(2+) 


7239 
7246 
7422 


7206 
7244 
7422 












{be) {us) 


^(2-^) 


7572 


7496 













system is significantly larger than in the corresponding heayy diquark (QQ)-light antidi- 
quark {qq) one. This fact is well known from the meson spectroscopy, where heavy quarkonia 
QQ are more tightly bound than heavy-light mesons Qq. For instance, we found that some 
of the {cu){cu) tetraquarks lie below open charm thresholds while all ground-state {cc){ud) 
tetraquarks are found to be above such thresholds. 

In Table |V] we confront our results for masses of heavy {QQ'){qq') tetraquarks with 



other theoretical predictions [2l|, [22, l23|, M, l25|, l26|, l27| . In Ref. [21|] the authors solve the 



four-body problem using the expansion in the harmonic-oscillator basis in the framework 
of the nonrelativistic quark model with a phenomenological potential. The same model 
with a different expansion basis, which can accommodate asymptotic states of two free 
mesons, is applied for the calculation of heavy-tetraquark properties in Ref. [13]. The 
stability of tetraquarks with heavy flavours is studied by using a variational approach and 

tetraquarks are analyzed 



a nonrelativistic potential model in Ref. j23|. In Refs. [2^ 
in the chiral constituent quark model using a variational approach. The potential of this 
model includes one-gluon, confinement and meson-exchange interactions. The existence of 
a virtual tetraquark state ccud is discussed in Ref. [26[ on the basis of semi-empirical mass 
relations. QCD sum rules are applied for the QQud tetraquarks in Ref. [23]. The main 
difference between our approach and the above quoted papers consists in that, from the 
very beginning, we explicitly reduce the relativistic four-body problem to the subsequent 
solution of two relativistic two-body problems assuming the diquark-antidiquark structure 
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of the {QQ'){qq') tetraquarks. From Table IVl we see that most of the presented approaches 
predict that only the l"*" state of the (bb) {ud) tetraquark lies below the open-bottom threshold 
(see also (Isl). In Ref. [22[ it is claimed that also the l"*" state of the {cc){ud) tetraquark is 
weakly bound against the DD* threshold, if it has a molecular structure. Note that such 
structures are absent in our approach. A large binding energy in the 1^ state of the {cc){ud) 
tetraquark is found only in Ref. [13] and it is claimed to originate from the meson-exchange 
part of the quark interaction potential. The recent QCD sum rule analysis [23] finds that 
only the {bb)(ud) tetraquark is expected to be a narrow state. 



V. CONCLUSIONS 

In this paper we have calculated the masses of the ground states of tetraquarks with 
two heavy quarks assuming the diquark-antidiquark structure. Such approximation allowed 
us to reduce the very complicated relativistic four-body problem to the solution of two — 
significantly more simple — relativistic two-body problems. All considerations were done 
in the framework of the relativistic quark model which proved to be successful in describing 
numerous properties of mesons and baryons [l, SB [o, 14, 15, 16, 17]. The parameters of the 



model were fixed previously from the meson sector and are kept unchanged in the present 
analysis. The diquarks were treated dynamically. Their masses and form factors, which 
take into account the diquark structure, were calculated on the basis of a numerical solution 
of the corresponding relativistic quasipotential equation. Note that they are the same as 
in our previous studies of light and heavy diquarks in heavy [sl and doubly-heavy baryons 
0], respectively. Light quarks and diquarks were treated completely relativistically without 
applying the v/c expansion. It was found that both the relativistic dynamics of light diquarks 
as well as their internal structure play an important role in the description of diquark- 
antidiquark bound states. The binding of a heavy diquark and a light antidiquark turned 
out to be weaker than the binding of a corresponding heavy-light diquark and heavy-light 
antidiquark. Thus, in contrast to the (eg) (eg) tetraquarks, which were discussed previously 
jsl, all the (ee)(gg') tetraquarks are predicted to be above the decay threshold into the open 
charm mesons. Only the I{J^) = 0(1+) state of {bb){ud) was found to lie below the BB* 
threshold. As a result, this state can decay only weakly and thus it should be narrow. The 
strange partner of this state {bb){us) is predicted to lie slightly (13 MeV) above the B*Bs 
threshold and, in principle, could be observed as a not too broad resonance decaying mainly 
into B*Bs. 

The investigation of the decay widths of heavy tetraquarks which are predicted to lie 
below the threshold of the open flavours represents another very important and interesting 
problem. It could be considered by means of the relativistic generalization of the analysis 
performed in Ref. 2^ . However, this problem is beyond the scope of the present paper and 



will be considered elsewhere. 
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